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Symmetric couplings among aggregates of n chromophores increase the transfer rate of excitons by a factor 
n 2 , a quantum mechanical phenomenon called "supertransfer." In this work we demonstrate how supertrans- 
fer effects induced by geometrical symmetries can enhance the exciton diffusion length by a factor n along 
cylindrically symmetric structures, consisting of arrays of rings of chromophores, and along spiral arrays. We 
analyse both closed system dynamics and open quantum dynamics, modelled by combining a random bosonic 
bath with static disorder. In the closed system case, we use the symmetries of the system within a short-time 
approximation to obtain a closed analytical expression for the diffusion length that explicitly reveals the su- 
pertransfer contribution. When subject to disorder, we show that supertransfer can enhance excitonic diffusion 
lengths for small disorders and characterize the crossover from coherent to incoherent motion. Owing to the 
quasi-lD nature of the model, disorder ultimately localizes the excitons, diminishing but not destroying the 
effects of supertransfer. When dephasing effects are included, we study the scaling of diffusion with both time 
and number of chromophores and observe that the transition from a coherent, ballistic regime to an incoherent, 
random-walk regime occurs at the same point as the change from supertransfer to classical scaling. 



I. INTRODUCTION 

In 1954 Dicke introduced the phenomenon of super radi- 
ance, a quantum interference effect induced by symmetries 
of spin-boson interactions, in which many interacting atoms 
could collectively conspire to yield an enhanced relaxation 
rate fl]. Due to this fundamental non-classical cooperation, 
the probability of a single photon emission from a n identical 
atoms collectively interacting with vacuum fluctuations be- 
comes n times larger than incoherent individual spontaneous 
emission probabilities \2\. Superradiant radiative relaxation 
can also occur in molecular aggregates due to inherent coher- 
ent feature of Frenkel exciton dynamics. This phenomenon 
can be observed when a closely packed group of molecules 
interacting under certain symmetry can collectively donate an 
excitation with a rate which is much faster than each indi- 
vidual molecule. The adversarial effects of inhomogeneous 
broadening and exciton-phonon interactions on such coopera- 
tive relaxation in molecular systems have been studied in de- 
tail su. 

The same symmetry principles that underlie superradiance 
can give rise an analogous phenomenon known as cooperative 
energy transfer or supertransfer EQ. Generally, the exciton 
transfer rate can be calculated from the transition probability 
of an excitation jumping from one molecule to another us- 
ing Fluorescence Resonance Energy Transfer (FRET), based 
on a perturbation treatment of the dipole-dipole interaction 
between individual molecules. However, under strong and 
symmetrized interactions of a group of n molecules the ex- 
citation becomes highly delocalized, leading to a large (effec- 
tive) dipole moment associated with the entire group. The re- 
sulting enhanced oscillator strength can lead to supertransfer 



when similar molecular assemblies, with comparable effec- 
tive dipole moments, exist that can play the role of acceptors. 
Under such conditions the rate of exciton dynamics should 
be calculated from these effective large dipole-dipole interac- 
tions to describe the coherent donation and acceptance among 
such molecular aggregates, with up to n? enhancement over 
the single molecule to single molecule transfer rate, even in 
the far field |9|. A primary goal of this paper is to study 
the behaviour of superradiance, supertransfer, and other non- 
classical collective phenomena in presence of disorder and en- 
vironmental fluctuations similar to those natural conditions of 
photosynthetic light-harvesting complexes PHFl l8l [lOl . Re- 
cently, the existence and role of quantum coherence in the 
dynamics of excitation energy transfer in biological systems 
have been of significant interest both experimentally ifTTHTTIl 
and theoretically 1(1814311 . Such studies can provide novel 
concepts and techniques for a deeper understanding of nat- 
ural/engineering excitonic systems potentially leading to effi- 
cient and robust artificial light-harvesting [29l [30l . 

A major problem in design and fabrication of novel exci- 
tonic devices is the limited exciton diffusion length that could 
be of about 10 nm in disordered materials. This issue has 
lead to low efficiency and complicated device structures in 
organic photovoltaic cells ll32l l33ll . and it is a bottleneck in 
the performance of excitonic transistors [34] and organic light 
emitting diodes [35 1. A key open question is whether one can 
use quantum-mechanical supertransfer effects to enhance ex- 
citon diffusion length in such disordered systems. Recent ex- 
perimental investigations of nano-engineered biological sys- 
tems suggest that under laboratory conditions certain aspects 
of photosynthetic complexes can be emulated that could be 
potentially exploited for efficient energy transport P6l [37l . 
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In particular, excitonic diffusion lengths up to a micron have 
been recently reported in engineered arrays of LH2 com- 
plexes [36|. Self assembled ring structures containing fluo- 
rescent chromophores attached to tobacco mosaic virus coat 
monomers exhibit efficient exciton transport together with 
a broad spectrum light collection with over 90% efficiency 
||37l . It is of considerable interest to explore whether the en- 
ergy transport mechanism in such systems was facilitated in 
part due to symmetries in the arrangement of chromophores, 
which give could rise to the phenomenon of supertransfer 0- 

m. 

In this paper, we investigate excitonic transport in systems 
consisting of rings of chromophores stacked in cylindrical ar- 
rays, as a function of the number of chromophores per ring, 
the spacing between rings, and the strength of decoherence 
and disorder. We also investigate excitonic transport in dipole- 
coupled spiral structures. Such geometries are relevant not 
only to the experimentally investigated systems mentioned 
above, but to naturally occurring cylindrical arrays of chro- 
mophores such as the green sulphur bacterium chlorosome. 
We use the symmetries of the system to derive analytic solu- 
tions for the behavior of the closed system in the absence of 
environmental interactions, and perform simulations to cap- 
ture the dynamics of excitonic diffusion in the presence of 
environmentally-induced noise and disorder. Our results pro- 
vide clear evidence for the presence of supertransfer in the 
appropriate regimes and for the destruction of supertransfer in 
other regimes. 

Specifically, we study the effect of supertransfer on the dif- 
fusion length a of an initially delocalized exciton along a lin- 
ear chain of chromophoric rings and helical rods, mediated 
by dipolar couplings between chromophores. Supertransfer- 
induced enhancements in the hopping rate translate into a 
commensurate increase in the distance travelled by the ex- 
citons. The strength of the supertransfer effect depends on 
the number of chromophores per ring n - higher n yields 
higher supertransfer rates, and on the distances between rings 
- smaller inter-ring distances yields larger asymmetries and 
diminishes supertransfer. By changing the strength of the 
interaction with the environment, we investigate excitonic 
transport in both coherent/ballistic and incoherent/diffusive 
regimes. By changing the degree of disorder, we investigate 
how the degree of localization depends on n. The analysis 
is carried out analytically in the closed system case, and nu- 
merically in the presence of disorder and interactions with a 
bosonic bath. Our analytic and numeric solutions allow us 
to investigate scenarios where rings are packed closely in a 
cylinder, and where chromophores are arranged in spirals. 



H. THEORETICAL MODEL 

A. A simple example 

The concept of supertransfer can be easily reviewed in the 
case of two spin systems A and B, each containing tia and 
riB sites respectively, invariant under permutation symmetry 
|9l . In the presence of such symmetric couplings, the hopping 



rate from the symmetrized single-excitation state of A to the 
symmetrized single-excitation state of B is tiatib times the 
hopping rate of a localized excitation from one of the A to 
one of the B sites. (The effect can be even larger - up to 
r? A nB - for multiple excitation states 0. For simplicity, this 
paper will focus on single exciton states.) Supertransfer can be 
captured in a simple way by the following symmetric hopping 
Hamiltonian: 

tia riB nA,n B 

i=l j'=l i=l,j=l 

(1) 

Here, i labels the tia sites, and j the tib sites of A and B, 
respectively. Introducing the angular momentum operators 

j n A n A Ub 

j a = 2E ct *' j a = E ct ±' J B=E <r ±' ^ 

i i i 

the previous Hamiltonian can be expressed as 

H = -e A J 3 A - £bJ 3 b + liJpB + JaJb) (3) 

From the above expressions, the probability amplitude for 
a single symmetrized excitation state \4>a) over the system 
A, \4>a) — S"=i to hop to a corresponding sym- 
metrized state over system B is given by jy/nA~nB, in first 
order perturbation theory. The corresponding transition prob- 
ability ^tiatib is therefore txatib times the probability j 2 
for a localized exciton state |1) to hop from one of the A sites 
to any one of the B sites. This enhancement effect in the hop- 
ping rate is called supertransfer |[7H3. The key feature of su- 
pertransfer is coherence within the individual systems A and 
B respectively. Supertransfer between A and B can either be 
coherent or incoherent, depending on the spatial separation of 
systems and strength of the interaction with the environment. 

B. Exciton supertransfer in cylindrical geometries 

Here, we study supertransfer effects in cylindrical aggre- 
gates of sites interacting via a 1/r 3 potential. We apply our 
analysis to the transition dipoles of an array of chromophores 
to look at excitonic hopping through the array. We consider 
diffusion in the context of two types of geometrical set-ups. 
In the first configuration, the sites are arranged in rings, which 
are stacked in a cylinder. Each ring has a radius R, con- 
sists of n dipoles each. There are a total of N rings, ar- 
ranged co-axially in a cylinder with a separation D between 
adjacent rings. The position of the i'th site inside the j'th 
ring is given by = (Rcos(2iri/n), Rsin(2iri /n) , jD), for 
i = 1,2, ... ,n and j = 1,2, ... ,N. In the second geome- 
try, chromophores form helical rods, with their positions given 
by Ti = (Rcos(2iri / n) , Rsin(2m/n), di/n), where d is the 
pitch of the helix, R its radius, n the number of chromophores 
per turn, with N turns in total, and i = 1,2, ... , nN . 

The single exciton manifold approximation will be em- 
ployed, so that the system Hilbert space is % = C nN , spanned 
by a basis S — {|m)}^. 1 , where |m) denotes the state in 
which the m th molecule is excited. 
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The Hamiltonian Ht of our ring aggregates is of the form 



H 



]((E)+ Vr )\r)(r\ 



E 

rn>r 



Jmr(\r)(m\ + \m)(r\), (4) 



with all the interactions included. For simplicity, we ini- 
tially consider an interaction of the form J mn = J/rf lm , with 
fmn — \r m — r n \. In what follows, the normalization J = 1 
will be adopted, with time measured in units of 1/ J. This in- 
teraction is a simplified version of the dipolar interaction that 
has the same distance dependence as the conventional dipolar 
interaction, but that does not take into account the orientations 
of the dipoles. This simplified assumption will be relaxed be- 
low. 

The influence of the environment over the system will be 
considered by including both on-site energy disorder and on- 
site dephasing (Haken-Strobl model) which has recently been 
widely used for studying environment effects on energy trans- 
port in light-harvesting systems, e.g., see Refs |[T9ll22ll24l . In 
(E) represents the average molecular excitation energy, 
which can be dropped by shifting all the energies by (E). The 
static inhomogeneous offset rj r in the energy of the r'th site 
reflects the disorder caused by the surroundings. We take r\ r 
to be a Gaussian random variable with standard deviation E 
and probability distribution given by 



P(r)r) 



1 



exp i 



(5) 



where we ignore correlations between the offsets of each 
molecule. The strength of the disorder £ is measured in units 
of J. 

The effects of the bosonic bath surrounding the system are 
be modelled by an on-site dephasing model or Haken-Strobl 
model, given by a Lindblad superoperator of the form 

C deph p(t) = 7 E [ S nP(t)Sn - S n p(t)/2 - p{t)S n /2] (6) 

71=1 

with the sum running over all the sites, S n — \n) (n\ and 7 the 
dephasing rate. Thus, the total dynamics is given by 



dp(t) 
dt 



- [H, p(t)] + £deph (/>(*)) - {ffrecom, p(t)}, (7) 



where the effects of exciton recombination is captured by the 
final term. Here, H Iccom — nJ2 n \ n )( n \> w i*h V K me life- 
time of the exciton. 



C. Choice of initial state and diffusion length 

By its definition, supertransfer relies fundamentally on two 
different symmetries: the permutation symmetry present on 
the Hamiltonian, as well as the permutation symmetry and co- 
herence on the initial quantum states. Both the ring chain and 
helical rods set-ups resembles the structures of the fluorescent 
chromophores attached to tobacco mosaic virus monomers, 
which could self-assemble into stacks of rings or helices (37). 
Note that the invariance under permutations of the n dipoles 



inside a ring resembles the symmetric interaction in ([TJ for 
D > R, where strong supertransfer effects should take place. 
When the rings become closer, D < R, then the interaction 
between rings becomes less symmetric. The strength of super- 
transfer is diminished, but the effect persists, as we will show 
analytically below. 

As the initial quantum state in the equations of motion, we 
choose a delocalized quantum state over the ring of n sites, of 
the form p Dcloc = \<j>) (<j>\, with \<p) = -±= Y^Li K)> supported 
on the middle (N/2) ring of the A-ring chain. We consider an 
odd number N of rings, each with an odd number of sites 
n, with AT = 2T + 1, n = 2t + 1. The rings are labelled 
by r — —T, . . . , -1, 0, 1, . . . , T, with denoting the middle 
ring. (The use of odd numbers of rings and sites are simply 
for the convenience of labelling.) To quantify how much an 
initial excitation located at the middle ring diffuses towards 
the boundaries, we solve for p(t) in the master equation ^ 
and calculate the second moment er(t) of the probability dis- 
tribution Pi(t) for an exciton to be at the i th ring at time t: 



a{t) = D 



\ 



(8) 



r=-T 



with p r (t) = Sj=i Pjjtf) th e probability for the exciton to 
be present in the r th ring at time t, Pjj(t) — (j\p{t)\j) the 
exciton population of site j on ring r and D the distance be- 
tween two adjacent rings. The diffusion length a(t) measures 
how far the exciton diffuses along the system of rings. We be- 
gin by presenting our analytic model for the closed quantum 
system to show that supertransfer enhances the diffusion of a 
delocalized initial state by a factor n. This enhancement for 
the closed system (coherent evolution and ballistic transport) 
will be found to persist in the open system, at least for short 
times. 



III. CLOSED SYSTEM DYNAMICS: ANALYTICAL 
CALCULATION OF DIFFUSION 

We will proceed to calculate a(t) for a closed system using 
a short time approximation, that is t\\H\\ -C 1, with \ \H\ \ the 
operator norm of the Hamiltonian H (maximum eigenvalue of 
H). We assume long lifetime 1/k t, such that recombina- 
tion effects can be ignored. An important observation is that 
the Hamiltonian of the linear chain of rings is very close to a 
block circulant structure. A block circulant matrix B € BC^.n 
is of the from 



B = circ(bo t>i 



-»N-1 



) = 



/ b 
bN l 



V bx 



bi 
bo 



b N 
b N 



b / 



(9) 

where the n X n matrices bi are themselves circulant, that is 
bj = circ(6j ; o h,i ■ ■ ■ &i,n-i) [38 1 . In this context, b rep- 
resents the Hamiltonian for a single ring, while the matrices 
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FIG. 1: (color online) Exciton diffusion in a chain of stacked rings as a function of number of nodes per ring n, per unit time and distance D 
between adjacent rings. A comparison is made between the diffusion obtained numerically (blue circles), and calculated via the theoretical 
formulas < | 1 4| > and \15\ (red stars). The upper (lower) panel shows diffusion obtained for an initial delocalized (localized) state with support 
over the middle ring of the chain. When the exciton is initially delocalized and the rings are far apart from each other as compared to their 
radius (the far-field regime D/R = 10), the set of inter-ring couplings are symmetric under permutations of the chromophores, and therefore 
diffusion grows linearly with the number of chromophores per ring n, as expected from supertransfer arguments. Meanwhile, if the exciton 
is initially localized over a single chromophore on the middle ring of the chain, the diffusion grows with the classical ^Jn scaling. When the 
rings are close to each other as compared to their radius (the near-field regimes D/R = 0.1 and D/R = 1), the diffusion has an anomalous 
scaling with the number of sites n, i.e. a ~ n a , a > 1, for both types of initial states. This is a finite size effect, as explain in the main text. 
We employed N = 31 rings. In all the cases, there is a good agreement between theory and numerical simulations. Note that for D/R = 0.1, 
both the time and coupling strength were adjusted in order to avoid boundary effects. 



t>i, i represents the n x n set of interactions between two 
rings. The circularity of bo is ensured by the symmetry of the 
ring structure. A symmetric, block circulant matrix B has the 
property that bj = bw-i, i = 0, 1, . . . , N — 1. This property 
would hold exactly if we were to impose periodic boundary 
conditions along the chain, so that the system is topologically 
a 'ring of rings' or torus. The actual chain Hamiltonian H, 
does not possess periodic boundary conditions. Nevertheless, 
for times short enough so that boundary effects can be ne- 
glected, we will see that the analytical expression for diffusion 
(7 obtained assuming block-circularity matches quite well the 
numerical results. 

A matrix B 6 BC^^ n if and only if it commutes with the 
unitary matrix II at „ = -K^®\ n , where 7Tjv is the fundamental 
N x N permutation matrix 7Tjv = circ(0 1 ... 0). 

Let us then consider the matrix H £ BCn,u as being 
block circulant, with circulant blocks. Such Hamiltonian 
will have the form H = circ(H Hi . . . Hx Hx • ■ • Hi), 
where Hj = circ(/ij j fe)^ =1 . The spectrum of such ma- 
trices is well-known |39|. Let Wj = exp(27rzj /N), and 
Pk = cxp(27ri/c/n). Define the matrix Fn(lj) given by 



1 W\ 



1 

K) 2 



(10) 



\1 W N -1 (w N -l) 2 ■■■ (W N -!) N 1 ) 

Then, the circulant matrix H can be diagonalized by Fn (w) ( 



F n (p), with eigenvalues e(p, q) given by 



E^Lo 1 YJk=l( w j) p ( w k) q hj,k 



(11) 
(12) 
(13) 



that is, it is the double discrete Fourier transform of the co- 
efficients hj t k, where the index j labels a block or ring, 
and k the element inside the j th block. For H real and 
symmetric, its real eigenvalues e(p, q) become T doubly 
degenerate for p = 1,2,...,T, and t doubly degener- 
ate for q = 1,2, ...,t, and are of the form e(p,q) = 
J2j=o Efc=o h Ptq cos(2-Kjp/N)cos(2irkq/n). The corre- 
sponding complex conjugated eigenvectors can be combined 
into two real eigenvectors. With this representation of the 
eigenvectors, the expression for the diffusion a for a finite 
chain of rings, for an initial delocalized state is given by 



CDeloc(i) = Dt^ 



n-1 



\ 3=1 



fe=0 



where D is the spacing between the rings. For a localized 
state, the corresponding expression for the diffusion is 



oloc(£) = Dt. 



T n-1 



\ 0=1 



(15) 



fe=0 
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FIG. 2: (color online) Exciton diffusion along a spiral chain of chromophores as a function of number of nodes per ring n, per unit time and 
pitch or aperture d. The spiral structure is approximated by a stacked ring geometry, in which each turn of the helix is replaced by a ring, with 
as many chromophores per ring as nodes per turn in the helix, and distance D between rings equal to the pitch d of the helices. The numerical 
results of diffusion obtained from solving eq.jTJ (blue circles) for the helical rod is compared with the diffusion along the corresponding 
approximating circular structure, using formulas l |14[ > and {15} (red stars). The upper panel shows diffusion obtained for an initial delocalized 
state over a set n of contiguous chromophores along a full turn of a helix, while in the lower panel an exciton with support over a single 
chromophore was employed. When the exciton is initially delocalized and the pitch of the helix is bigger than the radius (the far-field regime 
d/R = 10), the diffusion grows linearly with the number of chromophores per ring n as expected from supertransfer arguments. Meanwhile, 
if the exciton is initially localized over a single chromophore on the middle ring of the chain, the diffusion grows with the classical \fn scaling. 
When the rings are close to each other as compared to their radius (the near-field regimes d/R = 0.1 and d/R — 1), the diffusion has an 
anomalous scaling with the number of sites n, i.e. a ~ n a , a > 1, for both initial states. This is a finite size effect, as explain in the main 
text. The approximation of the spiral structure via a set of facing rings works very well, specially for the delocalized case. We employed 
TV = 31 turns/rings in our numerical simulation. Note that for D/R — 0.1, both the time and coupling strength were adjusted in order to 
avoid boundary effects. The values of diffusion are given per unit time and pitch d. 



These expressions are valid for small times or interaction 
strength, that is, t\\H\\ <C 1. The proof is given in the ap- 
pendix, for any initial condition. A couple of simple cases 
will help us understand these expressions. 

In the "far-field" regime, with D > R, the couplings among 
the rings will be highly symmetric, similar to ([T]i. For this 
case, t>i = bN-i = VJ n , with the rest of the block matri- 
ces t>i = 0, where J„ = circ(l 1 ... 1) is an n X n matrix 
with ones as entries and V = J(l/D nm ) 3 is the strength of 
the symmetric couplings between two nearest neighbouring 
rings, separated by a distance D. The expressions above for 
the diffusion will read (JDeioc = V^VDtn — \/2Jtn/ D 2 , 
&Loc — V^Jt^/n/D 2 , so that the symmetric couplings be- 
tween the rings produce a linear scaling of the diffusion as 
a function of the number of sites n, and yfn, for a delocal- 
ized and localized initial states, respectively. Thus, we proved 
that a linear dependence with the number of particles n occurs 
under symmetric couplings, and will be use in the rest of this 
work as a proof of supertransfer. We also obtain that for a 1/ r 3 
dipole-dipole decaying interaction, the diffusion decreases as 
1/.D 2 , with D the distance between rings, and is directly pro- 
portional to the typical interaction strength J. If all the inter- 
actions among the rings were included, the above expressions 
change to a De i oc = ^\f2Jtn/D 2 , a Loc = fg \f2Jt^/n/D 2 . 

In the "near-field" regime, with D < R, the scaling with 



n can increase beyond n or y/n for aoeioc and aLoc, re- 
spectively. The reason comes from how the interaction terms 

J2k=o hj.k and YskZo(hj,k) 2 scale up with n. Without loss 
of generality, let's consider two nearest neighbours rings. 
The interaction term for the delocalized case is of the form 
E£o%,* = E£o {D 2 + 2i? 2 (l - cos(27rfc/n)))- 3/2 , 
which exhibits a scaling with n faster than linear from D < R, 
for fixed j. This is a small size effect, and the scaling is linear 
for n large enough (as can be readily checked by approximat- 
ing the above sum by an integral for n ^S> 1). 



The expression for the diffusion of a delocalized state ( 14 1 
shows also some interesting "interference effects", in the 
sense that partial cancellations can occur when some of the 
coefficients hj k are negative. In particular, when the condi- 
tion 2fc=o — is met, the delocalized state does not 
propagate, while the localized does. Therefore, even in the 
closed system, an initial delocalized state does not necessarily 
diffuse longer than a localized state, during times t\\H\\ -C 1. 

Figure ([Tji makes a comparison between the simulated value 
of cr(t)/Dt, for delocalized and localized initial states, and 
for various values of the distance D between the rings. In the 
far-field regime D > i?, we clearly see a linear scaling of dif- 
fusion with n for a delocalized initial state, and with y/n for 
a localized state. For the near-field regime D < R, there is a 
superlinear scaling which arises from the behaviour of the in- 
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FIG. 3: (color online) Diffusion length a, along a system of stacked 
rings arranged co-axially in a cylinder, as a function of the number 
of chromophores n per ring and on-site energy disorder E that is 
represented in a logarithmic scale in units of J = 1. The system 
parameters are (iV, 7?, D) — (31, 1, 10). Each point corresponds to 
500 realizations. The time was set at t = 1. For small times or dis- 
order, tT, <C 1, we observe a linear scaling of diffusion with respect 
to the number of chromophores n, a clear signature of supertransfer, 
as the system evolves coherently. However, for long times or large 
on-site energy disorder, fE 3> 1, diffusive behaviour sets in and the 
mismatch among the on-site energies causes localization, regardless 
of the number of nodes per ring n. 



teraction terms in the Hamiltonian, as mentioned above. In all 
cases, there is a good agreement between the theoretical pre- 
dictions and the numerical simulations. For times longer than 
e/| | if ||, the excitons travel long enough so that the effects of 
the boundaries "kick in", and the block-circularity assumption 



Note that the Hamiltonian corresponding to the ring where 
the exciton is initially created does not enter equations ( 14 1 
and ( 15 1, as j = 0. Therefore, under the approximation 
< e, these expressions for diffusion do not depend on 



t\\H\ 

the couplings of the ring where the initial state is located. 

The above analysis, that was developed for a stacked ring 
geometry, can be used to approximate the diffusion in the he- 
lical rod setting, if the helices in the rod are approximated by 
stacked rings of the same radius R, with a distance between 
them equal to the pitch of the helix d = D, and with number 
of chromophores per ring n, equal to that in the rods per turn. 
Figure (|2]i exhibits this approximation. The blue circles cor- 
respond to numerical simulations for the spiral geometry with 
given R and pitch d, while the red stars depicts the formu- 
las <\\ty , { 15 i obtained from a set of facing rings with radius 
R and separation D — d, indicating that the approximation 
given by the rings structure is quite good. In principle, if the 
number of chromophores per turn is not an integer, the angu- 
lar position of the rings should be shifted by a given amount 
to account for this. However, due to the parametrization of the 
spiral that is being used, the nodes are perfectly aligned along 
the axis of the spiral, so no "shift" of the nodes position inside 
the rings is necessary. 




Disorder 



FIG. 4: (color online) Scaling exponent A of diffusion with time, 
a(t) oc t A , along a system of stacked rings arranged co-axially in 
a cylinder, as a function of the number of chromophores n per ring 
and on-site energy disorder E. The disorder is represented in a log- 
arithmic scale in units of J = 1 and the time was set at t — 1. The 
parameters are (N,R,D) = (31,1,10). The exponent A signals 
the character of the exciton propagation, with A = 1 representing 
quantum coherent or ballistic spreading, while A = 1/2 signaling 
classical or diffusive evolution. We observe that the exciton motion 
is coherent for small values of disorder, and shows a crossover from 
coherent to incoherent around E ~ 10 -2 , independently of the num- 
ber of nodes per ring n. For disorder E> 1/f the diffusion shows 
fluctuations around the classical random walk value A = 1/2. 



IV. OPEN SYSTEM DYNAMICS: NUMERICAL 
SIMULATIONS 

Having elucidated the behaviour of diffusion length in a 
closed system dynamics, in this section we analyse the be- 
haviour of supertransfer subjected to an open quantum system 
dynamics, via two effects: on-site energy disorder and on-site 
dephasing. 



A. Effects of on-site disorder 

For small values of static energy disorder, the motion of the 
exciton should be hardly affected, propagating coherently. For 
times t ^> X -1 , diffusive behaviour sets in and the mismatch 
among the onsite energies causes localization [40|. This is de- 
picted in Figure (|3j. For disorder £ l/t, the motion is co- 
herent, and we obtain a supertransfer regime a ~ n, while the 
exciton gets localized for ffi 3> 1. A similar conclusion can 
be drawn from Figure (j4j), where the exponent A in a(t) ~ t x 
shows a coherent (A = 1) evolution for small values of dis- 
order, and a crossover to classical diffusion takes place when 
tYi ~ 1/, independently of n. In the classical regime, the 
diffusion fluctuates around A = 0.5. 

Figure |5]) represents a level curve for diffusion as a func- 
tion of n and disorder S, obtained from Figure Q. For disor- 
der £ < 1CP 2 , the system is resilient to changes of disorder, 
and we see a linear enhancement of diffusion with n, while 
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FIG. 5: (color online) Level curve for diffusion a along a chain 
of stacked rings, as a function of disorder E and number of chro- 
mophores per ring n. This figure was obtained from Fig. {3}. For 
on-site energy disorder £ < 1CP 2 , the exciton still evolves coher- 
ently and its diffusion is independent of disorder, with a supertrans- 
fer scaling linear in n. Around E ~ 10 -2 we have a crossover from 
ballistic to diffusive behaviour, and a yfn type of behaviour sets in. 
This crossover agrees with Fig. Q. In the classical regine, for a 
fixed value of diffusion, the system is more resilient to disorder as n 
increases. 



for E > 10~ 2 , diffusion has a y/n behaviour. The areas be- 
tween different values of disorder provides an indication of 
how robust diffusion is under changes of disorder. From this 
observation, it derives that diffusion seems more robust un- 
der changes of disorder in the quantum or ballistic regime (for 
E < 10~ 2 ) than in the ballistic regime. 



B. Effects of on-site dephasing 

On-site dephasing diminishes the coherence among the 
sites of the i ring, which is fundamental for supertransfer ef- 
fects to take place. For a tight-biding Hamiltonian on an in- 
finite chain with periodic boundary conditions, nearest neigh- 
bour interaction J, and dephasing strength 7, see PP . the 
diffusion length is given by: 



— It -(1 

fi 2 7 L 7 



-ft\ 



(16) 



For small times, ft <C 1, the transport is initially ballistic, 
cr[t) ~ t, while for jt 1, the transport is diffusive, a(t) ~ 
\/t, with a crossover time of the order of 1/7. 

A similar behaviour should be expected for the diffusion for 
the ring structure. For small times compared to this dephasing 
scale, t < 1 /-f, supertransfer effects will enhance the diffu- 
sion length by a factor of n. This is shown in Figure Q, for 




FIG. 6: (color online) Diffusion length along a chain of stacked rings, 
as a function of number of nodes per ring n and time, in the presence 
of on-site dephasing. The dephasing strength 7 was set at 7 = 1, 
An initial delocalized state was employed for the numerical simu- 
lations. We observe that for times short compared to the dephasing 
time, t < I/7, the motion is still coherent and the diffusion is pro- 
portional to the number of chromophores per ring n, a(t, n) ~ n, 
a clear signal of supertransfer effects. Meanwhile, for long times 
or strong dephasing, t > 1 /y, the quantum coherence is destroyed, 
with diffusion exhibiting the classical scaling ait, n) ~ \fn. The 
parameters used are (N, R, D) = (31, 1, 10). The axis for a was 
multiplied by a factor of 100. 



7 = 1, where at small times the motion is still coherent and 
a(t,n) ~ n, while for t > 1/7 the motion is classical and 
a(t, n) r~j y/n. 

Using the block-circulant approximation employed in the 
closed system dynamics from the previous section, a similar 
calculation for diffusion can be carried out in the presence 
of dephasing as well. It is possible to verify that dephasing 
does not contribute to the diffusion, up to second order in the 
expansion of p(t) = e tc p(0), with C = — iadn + £d C ph, 
for p(0) being either a localized or delocalized state. That 
is, up to second order in time, the expressions ([T4| and ( p"5j ) 
remain unmodified under dephasing. Figure d7|) shows the 
diffusion as a function of the number of nodes n, for an ini- 
tial delocalized state (left panels), localized state (right pan- 
els), in the near field D/R = 0.1, D/R = 1 and far field 
D/R = 10. Most importantly, the theoretical predictions 
( pi] ) and 



15 1 (blue squares), are contrasted with numerical 
simulations of diffusion including many values of dephasing: 
7 = 0.1 (red triangles), 7 = 1 (green diamonds) and 7 = 10 
(brown squares). In all cases, for small values of dephasing 
7 = 0.1, the agreement between the numerical and analytical 
results is quite good, and the relative error between them is at 
most 3% for all cases. 

Figure ^ shows the best fit for the exponent A in a(t) ~ 
t x , for a fixed value of dephasing and different number of sites 
n. The exciton start initially with a ballistic spreading A = 1, 
decaying to a diffusive classical walk regime A = 0.5. As 
the number of sites increases, we see a sharper drop of the 
exponent towards 1/2. The dependence with n comes from 
the fact that the decoherence rate of a symmetric state coupled 
symmetrically to an common on-site dephasing bath goes as 
n times the single exciton rate 7. We lack an explanation for 
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FIG. 7: (color online) Effects of dephasing on diffusion. Exciton diffusion in a chain of stacked rings as a function of number of nodes per 
ring n, per unit time and distance D between adjacent rings. We compare the theoretical predictions for diffusion l |14[ > and < | 1 5[ > (blue squares) 
with numerical simulations under different values of on-site dephasing: 7 = 0.1 (red triangles), 7 = 1 (green diamonds) and 7 = 10 (brown 
squares). We employed TV = 31 rings and time was fixed at t = 1. The upper (lower) panel shows diffusion obtained for an initial delocalized 
(localized) state with support over the middle ring of the chain. When dephasing is included, it does not modify the theoretical predictions eq. 



14 1 and 



15 1, up to second order in the expasion of e p(0) in time. Therefore, for values of dephasing small compared with time (7 = 0.1) 



the agreement between the numerics and analytics is quite good. 



Scaling cr~f 




time 



FIG. 8: (color online) Scaling exponent A of diffusion with time, 
ait) oc i A , for an initial delocalized exciton moving along a system 
of stacked rings, as a function of time, for n = 1 (blue circles), 
to n — 6 (pink diamonds). On-site dephasing was included, with 
parameter 7 = 5. The exponent A signals the character of the exciton 
propagation, with A = 1 representing quantum coherent or ballistic 
spreading, while A = 1/2 signaling classical or diffusive evolution. 
We observe the exciton exhibits a ballistic spreading a{t) oc t (A = 
1) at small times ty <C 1, converging towards a classical random 
walk behaviour a(t) oc \[t (A = 1 /2) as t~y 2> 1, with the crossover 
around t ~ I/7. We lack an explanation for the dip of the scaling 
exponent A below 0.5 around t = 1, which may point to an intriguing 
physical feature, and is left for future research. The ring parameters 
were set to (TV, R, D) = (31, 1, 10). 
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FIG. 9: (color online) Scaling exponent a of diffusion with num- 
ber of chromophores n, cr(n) oc n a , for an initial delocalized exci- 
ton moving along a system of stacked rings, as a function of time, 
for on-site dephasing strength 7=1 (blue circles) up to 7 = 11 
(open brown circles). The exponent a provides an indication of su- 
pertransfer efffects when a = 1, while a — 1/2 signals classical or 
diffusive evolution. We observe the exciton exhibits a supertansfer 
regime cr(n) oc n (a — 1) at small times t*y <§C 1, converging to- 
wards a classical behaviour u(n) oc y/ri (a = 1/2) as t-y ^> 1. The 
crossover from one regime to another occurs around t ~ 1 /j, coin- 
cident with the transition from ballistic to diffusive regime in Figure 
{jj). The ring parameters were set to (TV, R, D) — (31, 1, 10), and 
we employed n = 1, 2, . . . , 10 nodes per ring in order to perform 
the fitting <r(n) ~ n a . 



the dip of the scaling exponent A below 0.5 around t = 1, 
which might point to an intriguing physical feature, and is left 
for future research. 

The scaling a in a(t) ~ n a as a function of time, for dif- 



ferent values of dephasing is shown in Figure ^9). The scaling 
goes from the supertransfer a = 1 to the classical walk regime 
a = 1/2 around the same time when the excitons diffusion 
starts to deviate from the ballistic spreading A = 1. 
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V. CONCLUSIONS 

We have analysed both analytically and numerically how 
symmetry -enhanced supertransfer, for circular and helical ge- 
ometry of chromophores, enhance excitonic diffusion lengths. 
For a closed system dynamics, we derived explicit expressions 
for diffusion for any initial state, demonstrating a factor n 
enhancement in diffusion when the inter ring couplings are 
similar among each other. These formulas explicitly capture 
quantum interference effects for initial delocalized states, in- 
dependently from the details of the rings where the excitons 
are initial created, and approximate the diffusion along helical 
rods quite well. 

Moreover, we have studied the effects of environmental in- 
teractions on supertransfer, by including both energy disorder 
and interactions with a bosonic bath in a Haken-Strobl pure- 
dephasing model. The n 2 enhancement in diffusion prevails 
for times smaller than the disorder strength, and the system 
becomes more resilient against random perturbations in the 
onsite energies for increasing values of the number of nodes 
n per ring or turn of the spiral. Due to the quasi ID nature of 
the aggregates considered here, supertransfer effects are frag- 
ile to disorder on the site basis energies. We have shown an- 
alytically that the dephasing does not affect the form of our 



closed expressions for diffusion, up to second order in time, 
in good agreement with our simulations. A numerical analy- 
sis enabled us to obtain the scaling exponents of sigma with 
number of nodes n and time t, revealing a crossover from the 
ballistic-supertransfer regime, to a diffusive-normal regime, 
in both cases around the same characteristic time 1/7. Our 
studies on exciton supertransfer dynamics presented here can 
be generalized for other complex quantum systems interacting 
with the non-Markovian and non-perturbative environments 
using the new techniques that have recently been developed in 
Refs. |29] ED. 
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Appendix A: Derivation of the closed expressions for the exciton diffusion length 



As shown in |TT[ ), the eigenvalues of a block circulant matrix are given by e(p,q) — EjLo 1 Y^k=o( w j) P { w k) q hj,k = 
J2j k exp(2iTrpj I N) exp(2iTrqk/n)hj ! k- Assume the initial state is given by "0(0) = J2^=i ET=o a r,s\r)\s), where for sim- 
plicity we assume a rs to be real. The second moment or diffusion a 2 (t) is given by a 2 (t) — D 2 E^=o R 2 Pn(t)> where pn{t) — 
Es=o I C^l {^\ l l ) {^))\ 2 i s tr, e probability that the exciton is in the ring R. Using the spectral decomposition of of the evolution op- 
eratatf(t)=ejq?(-ifff) we obtain A R , s {t) = (R\ (S\ifi(t)) = Ej^o exp(-ite(j, k))(R\e(j)) (S\e(k)) (e(j, k)\if)(0)), 

where we have used that the eigenvectors have the tensor product structure \e(j, k)) = \e(j)) \e(k)). For t\ \H\ | <C 1, we employ 
the approximation exp(— ite(j, k)) ~ 1 — ite(j, k). Putting all the terms together and expanding up to second order in time, we 
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, where we have used the convention that 



Es=o 8 Rt ol/n + 
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Therefore, the diffusion is given by <r(t) 2 = 



D 2 E«=-(T-i) ^ 



>R,0 



D2t2 J2 R =-(T-D r2 (J2 
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For a localized state on the 



ring and site 0, we have a Rt s = $raSs,o, and p R (t) ~ Es=o 5 R,o^/n + i 2 (Ej=o Y%=o h 3,k s -RJ s -S, 
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5R,ol-/n + t 2 (Y^=oYX=l h i^-RJ S -s,i 
D 2 Er=-(t-ii ^ 2 $R,o 



2 „_1 2 
1 l^S=0 n R,S 



= $R,0 
,T-1 



t 2 E 



h 2 

S=0 n R„S> 



and the diffusion is given by 



^ L Z^fl=-(T-1) n l~/S=0"'R,S- 
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